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ABSTRACT 

We present a theoretical analysis of the optical light curves (LCs) for short-period high-mass transiting extrasolar planet systems. Our 
method considers the primary transit, the secondary eclipse, and the overall phase shape of the LC between the occultations. Phase 
variations arise from (i) reflected and thermally emitted light by the planet, (ii) the ellipsoidal shape of the star due to the gravitational 
pull of the planet, and (iii) the Doppler shift of the stellar light as the star orbits the center of mass of the system. Our full model 
of the out-of-eclipse variations contains information about the planetary mass, orbital eccentricity, the orientation of periastron and 
the planet's albedo. For a range of hypothetical systems we demonstrate that the ellipsoidal variations (ii.) can be large enough to be 
distinguished from the remaining components and that this effect can be used to constrain the planet's mass. To detect the ellipsoidal 
variations, the LC requires a minimum precision of 10~ 4 , which coincides with the precision of the Kepler mission. As a test of 
our approach, we consider the Kepler LC of the transiting object HAT-P-7. We are able to estimate the mass of the companion, and 
confirm its planetary nature solely from the LC data. Future space missions, such as PLATO and the lames Webb Space Telescope 
with even higher photometric precision, will be able to reduce the errors in all parameters. Detailed modeling of any out-of-eclipse 
variations seen in new systems will be a useful diagnostic tool prior to the requisite ground based radial velocity follow-up. 
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1. Introduction 

Two observational methods have dominated the study of extraso- 
lar planets so far: radial velocity (RV) measurements and transit 
light curve (LC) analyses. Both have advantages and disadvan- 
tages. While RV determinations provide estimates of the plan- 
etary mass (M p ), the eccentricity (e) and the semi-major axis 
(a), they do not constrain the inclination (i) of the orbital plane 
with respect to the observer, thus only lower limits to M p can 
be determined. The transit method, on the other hand, provides 
information on i, the ratio of the planetary to the stellar radius 
(R p /R s ), and the duration of the transit (D). So far, only a com- 
bination of both strategies yields a full set of orbital and physical 
parameters for extrasolar planets. 

Currently, two space-missions are targeted at the detec- 
tion of transiting ex trasolar planets: CoRoT lau nched in 2006 
(iDeleuil e t al. 1997) and Kepler launched in 2009 dBorucki et alJ 
119971) . Their instruments are monitoring hundreds of thousands 
of stars, discovering thousands of planet candidates, whose RV 
follow-up will take many years (and may never be completed for 
the faintest candidates). Without RV follow-up, the most funda- 
mental parameter of an extrasolar planet, its mass, remains unde- 
termined from CoRoT and Kepler observations alone. The mass 
is the crucial parameter classifying an object as a planet, brown 
dwarf or a star. 

High-accuracy photometry has akeady been used for a num- 
ber of systems to show that the planetary thermal emission, as 
well as the reflection of the stellar light from the planet, are de- 
tectable. In particular, Welsh et al. (2010) report the discovery 
of ellipsoidal variations in the Kepler LC of HAT-P-7. This is an 



effect more commonly known from close stellar systems, where 
phase-dependent light variation arises from the gravitationally 
distorted stars. In HAT-P-7, the planet is close enough and mas- 
sive enough to induce the same effect. 

Other authors (lLoeb & Gaudill200llFaigler & Mazehll20Tlh 
have akeady presented treatments of reflected light, thermal 
emission, ellipsoidal variations and Doppler beaming. Our ap- 
proach to the problem includes two main differences: (i) we in- 
vestigate the sensitivity of our technique on the adopted phase 
function (Section 3.2), and (ii) we attempt to solve separately 
for the reflected light and thermal emission from the planet. For 
hot Jupiters around hot stars, the temperature of the planet in- 
creases dramatically and reflected light alone fails to explain the 
LC. 



In this paper, we investigate whether the presence of ellip- 
soidal variations in a planetary system can be used to place 
meaningful constraints on the mass of the planetary compan- 
ion. In Sect. |2]we describe the basic phenomena included in our 
simplified model of a planetary system's LC. In Sect.[3]we com- 
pare the magnitudes of the relative contributions and consider 
the likely range of planetary systems for which ellipsoidal vari- 
ations ought to be detectable with both CoRoT and Kepler. We 
reanalyze the HAT-P-7 LC and demonstrate that the ellipsoidal 
variations can be used to estimate the mass of the planet to within 
10% of the radial-velocity measured mass. We summarize our 
conclusions in Sect. [4] 
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2. Parametrization of photometric transits 



Standard models of LCs that have been used before the advent 
of space missions are based on a fla t curve out of transit and 
a limb darkening during the transit. Seager & Mall en-Ornelasl 
(2003) proved analytically that for each package of planet-star 
characteristics there is a unique LC. Analyses of high-accuracy 
data from space required a revision of this simple approach. For 
a high accuracy LC, the model should incorporate the reflected 
light from the planet, thermal emission, ellipsoidal variations 
and Doppler boosting, which deforms the overall shape of the 
LC, and the secondary eclipse (Fig.[TJ. 




Phase 



Fig. 1. The standard transit model includes the primary transit 
/i (dashed line centred at orbital phase z = 0), the secondary 
eclipse fi (dotted line around z = 0.5) and the reflected light 
from the planet / (solid line). The model is based on a HAT-7-b 
twin planet (Table|2]). 

In Fig.|2]we show the geometry of an arbitrary transiting sys- 
tem assuming an elliptical orbit. Let i be the angle between the 
observer's line of sight and the orbital plane normal, while the 
angle between the observer's line of sight projected onto the or- 
bit plane and the periastron is labeled (p. The star is in the center 
of the reference frame and d is the distance between the star and 
the planet; the distance between the star and the planet during the 
primary transit is denoted by d+ J p , during the secondary eclipse 
both bodies are separated by the distance d^L. 

To decode the geometry of the system from the LC we split 
it into three sub LCs. LC f\ describes the course of the primary 
transit, when the planet blocks the star's light, describes the 
form of the secondary eclipse, when the star blocks and plane- 
tary light, and fa the rest of the LC, when both star and planet 
contribute to the total amount of light. 

2.1. Transit and Eclipse 

The duration of the transit f\ is given by 



D =* 



PR* 



t/ PT 



i\ vr 



(1) 



where «d is the semi-major axis of the system, P is the or- 
bital period, R+ and R p are the radius of the star and the planet, 
respectively, e is the orbital eccentricity, b - d^i cos ///?* is the 
impact parameter dFordetaU 120081) and i is the inclination of 
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Fig. 2. This sketch of the transiting planetary system, as seen 
from above, explains the variables used in our calculations. We 
assume clockwise rotation. 



the orbital plane with respect to the observer's line of sight. To 
model the shape of f\ we use Eq. ([TJ and the limb darkening 
equation 



■j- = 1 - «i(l -fi) - m 2 (1 -/u) 2 , 



(2) 



with Mi and U2 th e linear and quadratic limb darke ning coef- 
ficients respectively (ICraretll2004t ISozzetti e t al. 2007), [i as the 
cosine of the angle between the surface normal and the observer, 
and Iq and 1^ as the intensities at the stellar disk center and at fi, 
respectively. Once the period is known from observations, one 
can fit the model to the observations to deduce R p , i, and 
d^ p . The transit of the secondary eclipse, fi, is fitted with the 
same model but without the effects of limb darkening. 



2.2. Phase-dependent light curve 

The shape of the LC between the transit and the eclipse (/) 
contains five contributions: (i) /+, a contribution from the stel- 
lar photospheric flux, which almost certainly varies with time, 
but which we assume to be constant throughout this analysis; 
(ii) / re f, a contribution from starlight reflected from the face of 
the planet, which is phase dependent; (iii) f± contributions from 
the intrinsic planetary thermal emission (split by day and night), 
which depends sensitively on a large number of factors, includ- 
ing the spectral range of the observation A, the irradiation tem- 
perature, and the atmospheric composition and structure; (iv) 
/ e ii, a contribution from ellipsoidal variations of the star itself 
due to tidal forces from the planet; (v) flux variations /d op due 
to a Doppler shift (Doppler boosting) of the stellar spectral en- 
ergy distribution with re spect to the bandpass of the instrument 
(Maz eh & Faiglerll20To|) . We normalize the LC by the stellar 
flux, which can be determined as the minimum flux observed 
at phase z = 0.5. The phase LC then becomes 



Mz) = 



f* + /ref + /th,d + /th,n + /ell + / 



dop 



/* 



(3) 



In the following, we study each of these terms in more detail. 
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2.2.1. Reflected light 



2.2.3. Doppler boosting 



The phase pattern of the reflected stellar flux depends on the 
phase angle z, i.e., the angle between star and observer as seen 
from the planet. Counting the orbital phase 9 from primary min- 
imum, the angles z, i and 9 are related through 



cos(z) = -sin(/)cos(0). 
The reflected flux / re f can then be expressed as 

.2 



(4) 



/ref(z)=a g /*(^)"0(z), 



where a g is the geometric albedo of the planet, /* is the stellar 
flux at a distance d from the star, and O(z) is the so-called phase 
function. It is not entirely clear which phase function provides 
the most appropriate description of extrasolar planets. A popular 
assumption is 



<J>l(z) = -(sinz + (n- z)cosz), 



(6) 



which models the planet as a Lambert sphere dRusselll[l916h . 
assuming that the intensity of the reflected light is = 1/tt at 
at phase z = n/2 (Venus case). An alternative choice is 



<J>c(z) = ^(l +cos(z)), 



(7) 



which assumes that the reflected light is <&c = 0.5 at phase z = 
n/2. In our model, we investigate both cases. 

The sum (p + 9 is related to the eccentric anomaly E through 



■ + 9 = 2 tan 



: tan (£72) 



(8) 



and E is related to the mean anomaly M through Kepler's equa- 
tion 

E = M-ecos(E). (9) 

If e + then the geometr ical albedo is not cons tant and it is 
given by the equation Eq. fTOldKane & Gelinoll2010l) 



„d-l 



A-d 



5(e d - 1 + e l - d ) 



3_ 

To 



(10) 



2.2.2. Ellipsoidal variations 

When a Jupiter-like planet orbits close to its host star, say 
a < 0.1 AU, then the star will be distorted. Hence, the sky- 
projected shape will vary, causing flux variations. Their magni- 
tude is given by 

/ell a M P (R^ ( l+ecos(<f> + 9) \ 3 . 3 . 



where 



/3 = 



\og l0 (GMJRl) 

logio^eff 



(12) 



is the gravity d arkening term (radiative energy transfer - 
Ivon Zeipell (1 19241) ). and M p and M* are the masses of the planet 
and the star, respectively. As these equations show, the detec- 
tion of the ellipsoidal distortion of the star provides information 
about the planetary mass with respect to the stellar mass. 



Finally, the relative flux because of stellar Doppler shifts given 

by 

^=(3-p)-, (13) 
with c as the speed of light, 

in on / sin(0) + e sin(0)\ 

K = 28.4 • P- 1/3 ■ M p • sin(z) • M* 2/3 — v ; (14) 

\ Vl -e 2 I 



(5) 

as the star's velocity amplitude, and p given by 



hcjkATax 



(3 - hc/kAT eff ) - 3 

? hc/kAT cS _ I ■ 



(15) 



Explanations of Eqs. ( fT3l - [151 ) and of the underlying phe- 
nomenon are given by Loeb & Gaudi d2003l) . If this Doppler 
shift can be measured in addition to the variations induced by the 
stellar ellipsoidal deformation, then indeed the planetary mass 
can be constrained independently from the stellar mass. 

In the next section, we will firstly explore constraints on ac- 
curacy raised by our method and then, as test cases, we will de- 
rive the mass of the known transiting exoplanet HAT-P-7, based 
on Kepler data. Finally, we predict the value of our method for 
the upcoming ESA mission PLATO. 

2.2.4. Thermal emission 

The equations for the planet ary emitted light are similar to those 
from lCowan & AgollOOlOl) . Thus 



and 



Si 



j e hc l AkT » - 1 

2 JiclAkT eS _ j 



s hc/AkT n _ I 



(16) 



(17) 



T e s is the effective temperature of the star, and 7a and T n are the 
temperatures on the day an d night side of th e pla net, respectively. 
Using the equations from |Hansen| (120081) and ICowan & Agol 
(1201 Oh . we calculate the mean temperature of the planetary day 
and night sides: 



Td = T Q ( l - or bo i) 



1/4 



12 



1/4 



T n = T (l- ff bo i) 



1/4 



(if- 



(18) 



(19) 



where T a = 7^(7?* /a!) 5 and < e < 1) is the energy 
circulation. 

The spectrum of the planet will strongly depend on the ab- 
sence or presence of clouds, the magnitude of the greenhouse 
effect and non-grey opacities. For the most of our analysis we 
assume that the planet emits as a black body. Anyhow, the dif- 
ferences between a black body and a non-black b ody emission 
are small (a factor of 1.5-3.5, see Hansen! 120081) . Assuming 
main sequence stars and planetary energy circulation e = 0.0 
(for the maximum thermal emission assuming black bodies), we 
found that for most hot Jupiters, the thermal emission from the 
planet is very weak in optical wavelengths. For our analysis we 
assume a blackbody, observed through the Kepler's sensitivity 
curve, which peaks at 575nm and covers 420 to 900nm (Kepler 
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Fig. 3. Thermal emission flux maps for three different bolometric 
albedo values (top : arboi = 0.0, middle a^oi = 0.3, bottom: aboi = 
0.5). Thermal emission is stronger for short term period planets 
around hot stars. For all the plots, we assume Rp = l.QRj. The 
white part of the diagram refers to /,/, < 1(T 7 . 



T e „ [K] 

Fig. 4. Color map of f± / f le f as a function of T e ff and Period 
for three different albedo cases (a g ,ab o i=0.05, a g ,arboi = 0.3, 
ag,a!boi = 0-5 for the top, middle and bottom figure respectively). 
We have fixed the energy circulation to e = 0.5 and Rp=l,0Rj. 
The white part of the diagram refers to fth/fref < 10 7 . 



Handbook). Figure [3] shows the star's effective temperature T e g 
versus the period of the planet for three different bolometric 
albedo (aboi) values. Colors refer to the expected stellar relative 
flux 

The thermal contribution only becomes comparable to, or 
even larger than, the reflected light component for short period, 



low albedo planets around hot stars. Distinguishing the relative 
contributions from a LC measured in a single passband, without 
any knowledge of the planet's albedo is hard, if not impossible. 
Both thermal and reflected components will have similar phase 
shapes. In Figure |4] we compare /,;, and f re f for three albedo 
cases (ar g ,aboi— 0.05, ttg,arb o i=0.3, <*g,a;boi = 0.5). 



4 



D. Mislis et. al.: Transiting exoplanet masses from precise optical photometry 



Table 1. Physical parameters of our two models. For the Super- 
Earth planet we assume no thermal emission (ffboi = 1 -0) 



Stellar parameter 


hot Jupiter 


Super-Earth 




1 CAD 

1.04«o 


0.8 1 K Q 


M, 


1.36 M e 


0.93 M e 




o/yu k 


52 ID K 


P 


0.53 


0.61 


6 


0.75 




Planetary parameter 


hot Jupiter 


Super-Earth 




1.16tfj 


0.15 Rj 


M p 


8.75 Mj 


0.0151 Mj 




0.30 


0.30 




U.UJ 


i n 

1 .u 


"l 


0.34 


0.20 




0.35 


0.57 


Orbital parameter 


hot Jupiter 


Super-Earth 


P 


5.63347 d 


0.85360 d 


i 


86.72° 


80.10° 


e 


0.52 


0.05 




72° 


5° 



The ellipsoidal variations occur on half the orbital period, 
and present a rather different phase dependent LC. Even if 
we cannot accurately distinguish between thermal and reflected 
components, we show in this paper, that they can be jointly dis- 
tinguished from ellipsoidal variations (Pig. [7J. 

3. Results 

3.1. Simulations 

To assess the data accuracy required by our method to yield sat- 
isfying constraints on the planetary parameters, we simulate two 
LCs using the equations from above. One of our model planetary 
syste ms is a hot Jupite r, analogue to the transiting planet HAT- 
P-7b dPal et al.l l2009bl) . whereas the other one rese mbles the 
transiting Super-Earth CoRoT-7b (lOueloz e t al. 2009). We cus- 
tomize these models in terms of the geometric al bedo, for which 
we op timistically apply a g = 0.3 in both cases (Sudarsky et al. 
2000). Though observations of CoRoT-7b are reconcilable with 
e — 0, we chose e = 0.05 to simulate enhanced phase variations 
in the LC. After all, we are not heading for a reconstruction of 
these systems but we want to estimate how accurate comparable 
systems could be parametrized and, as an example, if a putative 
small eccentricity of CoRoT-7b could be determined. This is par- 
ticularly important given the significant activity of CoRoT-7b. 

To these models, we add increasingly more noise simulat- 
ing data accuracies between 10~ 7 and 10~ 4 . The phase effect 
in the LCs, i.e. in f$, is significant only for accuracies < 10~ 4 , 
which is w hy this phenomeno n has not been detected in the LCs 
of CoRoT dCostes et alj|2004l) . We then fit the noiseless model 
from Sect. [2] to each of these - more or less - noisy LCs and 
use 1000 Monte Carlo simulations to calculate the standard de- 
viations for e ach parameter in eac h fit using a bootstrap method 
explained by Alon so et al.l d2008l) . In Fig. [5] we show the rela- 
tive errors resulting from these fits for the planetary geometrical 
albedo a g (cr a /a g ), eccentricity (cr e /e), orientation of perias- 
tron (o-^/ai) and the mass ratio (cr M /M) as a function of the root 
mean square (R MS) of the data. With an accuracy provided by 
Kepler of 10~ 4 dKoch et al .11201 Ot) , the eccentricity of a CoRoT- 
7b-like Super-Earth could merely be determined with a unsatis- 
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Phase 

Fig. 6. Residuals flux between Lambert & geometrical sphere. 
The residuals curve is identical to the ellipsoidal variations 
curve. By selecting the wrong phase function we over/under es- 
timate the planetary mass. 



factory relative error higher than 1.0 which has no physical use. 
For a hot Jupiter similar to HAT-P-2b, however, the relative er- 
ror in e is about 0.001 in the best case (hot Jupiter - RMS 10~ 7 ). 
The orientation of periastron for the Super-Earth could be con- 
strained to approximately +4° in the worse case (RMS 10 -4 ) 
while for the hot Jupiter the accuracy is as low as ±0.001° in the 
best case (RMS 10~ 7 ). Restrictions of the geometrical albedo a„ 
are +0.001 in the best case (hot Jupiter - RMS 10~ 7 ) to + 1.1 in 
the worst case (Earth like - RMS 10~ 4 ) with no physical mean- 
ing. The respective accuracies in planet mass Mp are ±0.028 Mj 
for hot Jupiters (best case) and ±0.95 M/for Super-Earths (worst 
case). 

3.2. Relative Contributions 

Different phase functions (eqs. (O & (|7])) lead to different peaks 
in the thermal emission and reflected LC. The ratio between 
thermal emission and reflected light over ellipsoidal variations 
becomes smaller (for the Lambert sphere) or larger (for the ge- 
ometrical sphere), which is why different phase functions lead 
to different mass values. The mass difference (AMp) between 
the two cases is a constant (36.4%). For example, if we use a 
wrong sphere, the mass we measure for HAT-P-7b exoplanet is 
AMp ~ 0.67 M j smaller (geometrical) or larger (Lambert) than 
the true mass. Seeking out the most plausible phase function, we 
fit two models, both differing only in the assumption of either 
the Lambertian or the geometrical sphere, to the data. For most 
systems, if the components f re f, f,i, & f e u are very weak, it is im- 
possible to distinguish between the two phase functions. In this 
case we select the Lambert sphere in order to calculate the upper 
mass limit of the planet. In Fig.[6]we show the residuals between 
the two different phase functions. Notice that the shape and the 
phase (x-axis) of residuals is identical to ellipsoidal variations 
curve. 

For most hot Jupiter systems thermal emission in optical 
wavelengths is much weaker than variations due to the varying 
projected shape of the star. For orbital periods of 1 and 5 days of 
a GOV star, ellipsoidal variations are 3 and 10 times larger than 
thermal emission, respectively, while / t h//* ~ 10~ 5 and 10~ 9 
(Fig. [3] & HJ. As the star becomes hotter (F or A type), or the 
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Fig. 5. Relative errors in <?, a>, a g and Mp as functions of data accuracy. The solid line denotes the HAT-P-2b twin while the dashed 
line labels the CoRoT-7b analog. 



period shorter (< 5 days), the thermal component become to be 
more significant. 

Thermal emission, stellar ellipsoidal deformations, and 
Doppler shift changes are very weak phenomena (fainter than 
1CT 5 ), which require very high-accuracy observations to be de- 
tected. So far, it has been impossible to detect these effects from 
the ground but based on data from the recent space-based mis- 
sion Kepler, our system of equations offers a new gate for exo- 
planet mass determinations. Neither RV measurements nor tran- 
sit timing variations are necessary for our procedure. 

Our parametrization of reflected and thermal planetary emis- 
sion are nearly identical, so our model is by definition unable to 
distinguish between these two effects. Both equations IBl and [T6l 
are only function of the orbital phase (<£>l.c) and both compo- 
nents (thermal and reflected) provide exactly the same LC. We 
are not able to distinguish between the two components by fit- 
ting our model in a hot Jupiter system, which the planet is tidal 
locked by its host star. If the planetary rotation and orbital pe- 
riod are different (for a case of a high eccentric orbit and the 
simple assumption that there are no winds), then reflected light 
and thermal emission LCs will show different shapes. 

3.3. Application to Kepler data: The Case of HAT-P-7b 

Kepler is a space-based mission, launched on Mar. 7, 2009, with 
the primary goal to detect Earth-sized and smaller planets in the 



habitable zone of solar-like stars. For a my = 12 star, the photo- 
metric precision is lo gi n (RMS) — '4-35 f or 6 hr combined dif- 
ferential photometry dBorucki et alJl20lTh . We choose the well- 
studied transiting planet HAT-P-7 in the Kepler field, with mass 
and system parameters constrained by both photometric and RV 
observations, to test our photometry-only method. We use all the 
short exposure data (1 minute exposures) from Kepler's public 
archive. Our sample includes 67 transits and the phase folded LC 
was made by ~ 180, 000 epochs binned by a factor of ~ 50. 

IWelsh et"aT1 (1201 Ol) already found ellipsoidal variations in 
HAT-P-7. Using the established physical parameters of the sys- 
tem, we calculate that the effect of (/ e n ~ 3.1 ■ 10~ 5 ) is 
~ 3.56 times more pronounced than the Doppler shift variations 
(/do P ~ 8.7- 10" 6 ). Nevertheless, we include both effects in our 
model. 

We assume that the stellar radius and effective temperature 
are known (r e ff = 63507T). We fix the ratio of stellar and plan- 
etary radius R p /R+ = 0.0778, the inclin ation ; = 83.1° an d the 
relative semi-major axis a/R^ = 8.22 dWelsh et al.ll2.010h . For 
our first fit, we use Eq. [3] but we exclude f e u and fd„ p . We only 
fit f re f plus frh in order to check if only these two components 
could successfully explain the data Cf^ f+th ). Our second fit in- 
cludes all four components from Eq. [3] (X^n)- The x 1 ratio be- 
tween the reflected light plus the thermal emission model and the 
model which includes all four components is xtf+^/x^n = 3-04. 
Using an F-test, we found that the complete model explains the 
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Fig. 7. LC of HAT-P-7. The solid line shows the best fit solution 
of our synthetic model, i.e. including stellar ellipsoidal distor- 
tions and Doppler effect variations. The dashed red line refers to 
the reflected light component, the blue dotted line refers to the 
ellipsoidal variations and the purple, almost consant dashed line 
refers to the thermal emission. The Lambert sphere assumption 
produced the best fit. 



Table 2. HAT-P-7b system parameters. 



HAT-P-7 parameters 

e 
d> 



Our model 

1.27 ± 0.10 Mil M Q 
0.0 ±0.1 
0.0 ± 1.0° 
0.21 ± 0.08 



Other studies 

1.20 + 0.05 MjI M q ' 
0.003 ± 0.012" 
0.0" 
0.18" 



a (Welsh etal. 2 010) 
" (Pal et al. 20~9a) 



data better than the simple model of the reflected plus thermal 
emission. There is a 0.01% chance that the more complex model 
do not produce a better fit to the data (F-Ratio: 1 16.28, P- Value: 
0.0001). In Fig.|7]we show the observed LC of HAT-P-7 with our 
best fit as well as a comparison to a model, which includes only 
reflected stellar light from the planet, and our enhanced model. 
In Table |2] we list the output parameters of the fit. Because /,/, is 
very weak (10~ 6 - 10~ 7 , see Fig.[3]for the black body case), ar D ob 
Td and T„ values can no t be calculated accurately (very large 
errors) . IWelsh et all (1201 Ol) . used the secondary eclipse depth of 
HAT-P-7b, in order to measure the flux from the night side of 
the planet (/ f / v , = 2.2 • 10~ 5 ). From Fig. [3] (bottom diagram) is 
clear that the planet is much brighter and much hotter than ex- 
pected, and other heating mechanisms might exist. We can not 
disentangle f re f from /,/,, however this does not affect our main 
result. 

We find the ratio of the planetary mass over the stellar mass 
to be M D /M * - 1.27Mj/M , co mpared to M p /M* = 1 .20 M,/M Q 
deduced by IWelsh et alJ d2010h . Because the f e u component is 
uncorrected with f re f and/or /,/,, the accuracy of the mass ratio 
(Mp/Mi,) can be calculated with high accuracy. Although out- 
values are not as accurate as those from RV measurements, we 
can infer the planetary nature of the companion. 



3.4. PLATO mission 

PLATO is a prospective ESA mission aiming at the character- 
ization of transiting exoplanets. The mission will use 40 small 
telescopes covering 42 x 42 de g 2 in total wit h a supposed over- 
all accuracy of ~ 2.8 • 10~ 5 dClaudill2010T) . This will satisfy 
the needs of our procedure. In the following, we demonstrate the 
capabilities of our photometry-only method in parameterizing 
hypothetical exoplanet systems. 

We investigate the detectability of f e \\ and /dop in a range 
of different systems, assuming that our target stars belong to 
the main sequence. For each combination we calculate the am- 
plitude of the flux variations AF in the LC. In Fig. [8] we 
present three plots, each belonging to a different planetary mass: 
M p = 0.5 Mj (top), 1 Mj (middle), and 5 Mj (bottom). For com- 
parison we indicate some known transiting systems with similar 
masses. In the upper panel AF is dominated by f e \\. Here, for the 
sample of known systems AF > 10 , which would be chal- 
lenging to be detected. Nevertheless, with PLATO one will be 
able to detect f e \\ for planets with P < 3 d around stars with 
spectral type K0 or later (blue zone in the top panel of Fig. [8). In 
the middle diagram (M p = lMj) both / e u and /dop affect the LC 
and AF ~ 5 ■ 10~ 5 for most of the known systems - detectable 
with PLATO. In the bottom panel, where the model planet is 
most massive, the LCs are mostly affected by the Doppler shift 
variations rather than by the ellipsoidal geometry of the star. The 
phenomenon is fairly detectable for the most cases of the known 
planets while AF ~ 12- 10~ 5 . 



4. Conclusions 

The mathematical tools presented in this article can be used for 
a complete parametrization of transiting exoplanet systems on 
the basis of high-accuracy LCs only. In our model, for a specific 
range of systems, RV measurements are not necessary to con- 
strain the mass of the planet (M p ), orbital eccentricity (e), the 
orientation of periastron (a>), and the geometric albedo of the 
planet (a g ). Our model also incorporates the characterization of 
the ratio of planetary and stellar radius (R p /R+), orbital period 
(P), and the orbital inclination (/). In order for this method to be 
fully applicable, the planet must be massive enough and orbit its 
host star close enough as to distort the stellar structure signifi- 
cantly. 

With the current Kepler and potential future missions such 
as PLATO, we are able to measure the mass of the hot Jupiters, 
where the orbital periods must be < 1.5, < 4, and < 10 days 
in order to measure planetary masses of 0.5 Mj, 1 Mj, and 5 Mj 
respectively. Furthermore, for host stars with spectral types ear- 
lier than FV, thermal emission component is detectable in the LC 
(assuming black bodies). This will affect the accuracy and even 
the detection of ellipsoidal variation. 

As we show, the Kepler mission provides an accuracy suit- 
able enough for our procedure to be applied for the characteri- 
zation of extrasolar planets. We confirm the planetary natures of 
the hot Jupiters HAT-P-7 from Kepler data alone. Our method 
yields masses of M p = 1.27+° '° M s x M*/M for HAT-P-7b. 
Our results ( M P , a g , o-boi, e, Td, T„) are in good agreement with 
IWelsh et"aTI(l20 10) and we also confirm the ellipsoidal variations 
of HAT-P-7b system. Using ellipsoidal variations we calculate 
the planetary mass from the photometric LC itself, without any 
RV measurments. Our technique will benefit from future space 
missions such as PL ATO and the James Webb Space Telescope 



(iDeming et al.l l2009) with RMS 
tively. 



3 ■ 10 3 and < 10" 6 , respec 
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Fig. 8. Combination of stellar ellipsoidal distortion and Doppler 
shift effects in simulated PLATO observations. The panels show 
the amplitude of flux variations AF in the LCs for three differ- 
ent planetary masses: Mp = 0.5 Mj (top), 1 Mj (middle), and 
5 Mj (bottom) around a range of main-sequence stars. Diamond 
symbols represent known transiting planets with similar masses. 
The white part of the diagram refers to f e u < 10~ 55 . 
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